Spin freezing and magnetic inhomogeneities in bilayer manganites 
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We have performed a muon spin rotation study on poly crystalline samples of electron-doped lay- 
ered manganites, La2-2a;Sri+2 a ;Mn207 (0.4 < x < 1), in order to investigate the local magnetic 
structure and spin dynamics. Our results provide evidence for phase separation into A-type antifer- 
romagnetic and charge-ordered phases for x—0.52 and spin freezing at low temperatures (T < 100 K) 
for 0.52 < x < 0.75. A new phase diagram which includes this spin freezing region is proposed. 



The phenomenon of phase separation is crucial for de- 
termining the intrinsic properties of many magnetic ox- 
ides 0. The origin of the colossal magnetoresistance 
(CMR) in manganese perovskite structures is often linked 
to the coexistence of ferromagnetic (FM) metallic and an- 
tiferromagnetic charge ordered (CO) insulating domains 
or to the formation of nanoscale magnetic clusters 0, • 
One technique which provides information on the local 
magnetic order is muon spin rotation/relaxation (/iSR). 
In underdoped cuprates, /xSR has identified the coexis- 
tence of superconducting and antiferromagnetic phases 
as well as spin freezing believed to be related to the gen- 
eration of stripes . 

In this letter we provide direct evidence for spin- 
freezing and phase separation in electron-doped layered 
manganite structures using /iSR. We investigate the lo- 
cal magnetic structure and spin dynamics across a se- 
ries of magnetic structures of Lag— 2a:Sri-|-2xMn207 with 
0.4 < x < 1. For compounds in the Mn 4+ -rich half of the 
phase diagram (0.5 < x < 1) neutron diffraction studies 
[5j reveal a progression of antiferromagnetic (AFM) in- 
sulating phases: from A-type (0.5 < x < 0.66) to C/C*- 
type (0.75 < x < 0.9) and to G-type (0.9 < x < 1.0). 
Furthermore, no long-range magnetic order (NLRO) has 
been observed between the A-type AFM and C/C*-type 
AFM regions (a;=0.66-0.74) [5j. Our experiments mea- 
sure the temperature dependence of the order param- 
eter in the samples with long-range order (LRO) and 
show how this is affected by the interplay between mag- 
netic and charge order. Moreover, by cooling a sample 
in the NLRO region, we are able to follow the progres- 
sive slowing down of magnetic fluctuations and demon- 
strate the development of short-range order. Based on 
our results we are able to propose a new phase diagram 
of La2-2ajSri + 2a;Mn207 (0.5 < x < 1), supplementing 
the information already obtained from neutron diffrac- 
tion 0. 

Zero-field /iSR is especially suited for study of short- 
range magnetic correlations since the positive muons 
are a sensitive local probe. Because the local magnetic 
fields at muon sites result primarily from dipolar interac- 
tions which decay very quickly with increasing distance 
(~ 1/r 3 ) the effective range investigated with muons is 



~ 20 A. In a system having magnetically ordered and dis- 
ordered regions the /xSR data is composed of two differ- 
ent signals corresponding to different environments, with 
signal amplitudes roughly proportional to their volume 
fraction. In the case of a magnetically ordered system 
our data are fitted using a relaxation function, G z (t), of 
the form 

G z (t) = A r u exp(-t/Ti) + A osc cos(27W M t) exp(-t/T 2 ) (1) 

where the first term is the longitudinal (spin-lattice) re- 
laxation, the second term is the muon precession at fre- 
quency v^^^B^jlix in the local internal field, B^, in the 
ordered state, X/T% is the transverse relaxation rate and 
7 M /27r=135.5 MHz/T is the gyromagnetic ratio. The lon- 
gitudinal relaxation rate is defined by 1 jT\ oc (7^ Ai? M ) 2 r 
for rapid fluctuations, where AB^ is the amplitude of the 
fluctuating field and r is the Mn-ion correlation time. 
Near a magnetic phase transition r (and hence 1/Ti) in- 
creases with cooling due to the critical slowing down. For 
polycrystalline samples, where the local magnetization 
points along the muon spin direction with probability 
1/3, the corresponding amplitudes are = G z (0)/3 
and A osc = 2G 2 (0)/3. Eq. Q can be generalized in a 
straightforward manner if there are multiple frequencies 
in the data corresponding to magnetically inequivalent 
sites. Polycrystalline samples (~ 0.5 g) were prepared as 
described in Ref. 0. Zero-field /LtSR data for 10 com- 
positions were taken on the EMU beamline at the ISIS 
facility (Rutherford Appleton Laboratory, UK) and for 5 
compositions on the GPS spectrometer at the PSI muon 
facility (Paul Scherrer Institute, Switzerland). 

Fig. H shows /iSR spectra collected at T=1.8 K (at 
PSI). Two oscillations can be distinguished for the G- 
type AFM sample, 27=0.98, and only one damped os- 
cillation for all others, although for £=0.68 the damp- 
ing is extremely large. The data were fitted using Eq. 
and the obtained parameters are listed in Table [I] The 
source of the internal field is attributed to the aligned 
manganese moments in the vicinity of the muon. As 
a function of hole doping, x, the values of v ^ indicate 
changes in the local field when the spin structure varies. 
The precession frequency for x=0.40 is 1^=28(1) MHz 
at T = 1.8 K (which corresponds to an internal field of 
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FIG. 1: Muon spin relaxation data at T=1.8 K for 
La 2 -2xSri +2a; Mn207 with (a) x=0A0 (FM), (b) x=0.52 (A- 
type AFM), (c) a;=0.68 (NLRO), (d) 2=0.82 (C/C*-type 
AFM) and (e) 2=0.98 (G-type). Two oscillation frequencies 
can be observed for 2=0.98 and a very strongly damped os- 
cillation for 2=0.68. The solid lines are fits to the data using 
Eqd 



TABLE I: Parameters obtained by fitting the tiSR spectra for 
La2-2xSri +2 xMn20 7 at 1.8 K (Fig.[TJ. 



Sample i/ M 


1/T 2 


1/Ti 


l/(27iVT 2 ) 


(MHz) 


(MS- 1 ) 


(MS" 1 ) 


(%) 


x=0.40 28(1) 


57(6) 


0.005(1) 


0.32(4) 


x=0.52 68(4) 


110(14) 


0.10(1) 


0.26(4) 


x=0.68 30(5) 


150(23) 


0.06(2) 


0.80(18) 


x=0.82 5.4(2) 


20(1) 


0.03(1) 


0.58(15) 


x=0.98 47(2)/76(2) 


16(2) 


0.065 


0.054/0.034(0.6) 



oscillations is dominated by a static distribution of in- 
ternal local fields determined by the neighbouring Mn 
spins at the muon site ||. Nevertheless, some dynamic 
broadening, i.e. fluctuations of Mn spins, is also present 
since 1/Ti is nonzero. The relative width of the field 
distribution can be parameterized by \ j ^T-i) (Ta- 
ble HJ which is very large (« 0.80) for the NLRO sample 
(x=0.68) and quite small (< 0.055) for the G-type AFM 
compound (x=0.98) that has only 2% Mn 3+ ions mixed 
with Mn 4+ ions. A large relative width can be associated 
with a large distribution of Mn 3+ and Mn 4+ ions (result- 
ing in a large degree of local magnetic inhomogeneity) , 
or with the coexistence of different magnetic phases. 
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FIG. 2: Temperature dependence of the zero-field muon pre- 
cession frequency, v^, and internal field, B^, for (a) 2=0.52 
and (b) 2=0.98. The inset shows the temperature dependence 
of the relaxing asymmetry, A r i x for the two compositions. 



£^=0.21(1) T), similar to that found for the ferromag- 
netic system, a;=0.3 Q- The observation of two muon 
frequencies for the G-type AFM compound, x=0.98, sug- 
gests that the muon occupies magnetically inequivalent 
sites in this compound. In oxides the muon site is usually 
about 1 A from an oxygen atom Q . In the layered man- 
ganite structure there are three inequivalent oxygen sites 
0(1), 0(2) and 0(3) and one likely muon site is close to 
0(2) (0, 0, 0.2) 0. 

The longitudinal relaxation rate, 1/Ti, at T = 1.8 K 
is 2-3 orders of magnitude smaller than the transverse 
relaxation rate, I/T2, indicating that the relaxation of 



The temperature dependence of f M (T), (and hence 
S M (T)) for £=0.52 and x=0.98 is shown in Fig. For 
cc=0.52 we find a discontinuity in f M (T) around 100 K. 
This arises because of the presence of two magnetic 
phases, A-type AFM and a CO phase similar to that 
found in x=0.5 (lcHllj . which compete with each other 
in the region between 100 K and 200 K. Below 100 K 
the internal field increases significantly due to an addi- 
tional magnetic transition. This transition coincides with 
a broad maximum in 1/Ti jl2j . signifying a freezing of 
the local spins. This effect is strikingly similar to that 
observed in a bilayer cuprate 4], in which a freezing of 
spins of doped holes is superimposed on AFM long-range 



3 



ordered Cu 2+ spins. The slow variation of the relax- 
ing asymmetry, A r i x (inset of Fig- EJ between G z (0) (at 
200 K) and G 2 (0)/3 (below 100 K) further confirms the 
two-phase model whose volume fraction is temperature 
dependent and that the whole sample is magnetic below 
100 K. The behavior of the x=0.52 compound is in con- 
trast to the G-type AFM compound, a;=0.98, for which 
the drop in the initial asymmetry is much sharper (see 
inset Fig. 2(b)) and both precession frequencies ^i(T) 
and 1^2 (T) can be fitted to v^{T) = z^(l - T/T N f 
(Fig. 2(b)). This procedure, although usually applicable 
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FIG. 3: Temperature dependence of different fitting param- 
eters for 2=0.68; (a) the zero-field cooled (open circles) and 
field-cooled (closed circles) magnetization in 10 mT; (b) the 
static amplitude (A osc for T < Tf and A s for T > Tf) and 
the dynamic amplitude (A Tlx for T < Tf and A t for T > Tf); 
inset shows f M (T) and Akt/(27t); (c) longitudinal relaxation 
rate, 1/Ti, (the inset shows the transverse relaxation rate, 
I/T2, and slow component, A s ). The freezing temperature 
Tf (vertical dotted line) separates the static and the dynamic 
regions. The vertical dashed line shows Too defined as the 
maximum in the magnetization data, in agreement with elec- 
tron diffraction studies [13. The solid lines are a guide to the 
eye. 

only to the asymptotic critical regime, gives a reason- 
able description of the data over the entire temperature 
range (see Fig. H{b)), yielding Tn=145.0(5) K in agree- 
ment with neutron powder diffraction @ and £=0.24(2) 



(close to the value corresponding to a 2D XY system 
(/?=0.23) 00). 

One of the most intriguing regions of the Mn + -rich 
phase diagram is x=0.66-0.74 where NLRO has been de- 
tected using neutron diffraction Data for x=0.68 at 
T=1.8 K shown in Fig. 1(c) were analyzed using a re- 
laxation function which was (1) Eq. 1 for T < 70 K, (2) 
G z {t) = A s e~ Xst + A f e~ Aft for 70 < T < 120 K and 
(3) G z (t) = Aie- Xit for T > 120 K and the results of 
the fitting are presented in Fig. [21 At low temperatures, 
the strongly damped oscillation frequency decreases upon 
warming from ~ 30 MHz at 2 K to ~ 10 MHz around 
65 K (inset to Fig.^b)), suggesting the development of 
a quasi-static field at the muon-site indicative of short- 
range order. In the low temperature regime we found 
that modelling data using a Gaussian Kubo-Toyabe re- 
laxation function parameterized by Akt/(27t) gave sim- 
ilar results (see inset of Fig. Etb)). 

The broad static distribution of local fields is quantified 
by the large 1 /T 2 values which is replaced by the slow dy- 
namic relaxation rate A s above 70 K (inset to Fig. 3(c)). 
The damped oscillations disappear above 70 K, the tem- 
perature at which 1/Ti exhibits a pronounced maximum, 
corresponding to the slowing-down of spin fluctuations, 
and we identify this temperature as the spin-freezing tem- 
perature, Tf. The amplitudes of the two components ^4 sc 
and A r i x shown in Fig. 3(b) provide information about 
the volume fractions of the sample where muons expe- 
rience static or dynamic spin-lattice relaxation, respec- 
tively. Above Tf the oscillatory component, A osc , (and 
the corresponding frequency v^) disappears and is re- 
placed by a slowly fluctuating fraction with amplitude A s 
(Fig. Elb)). This fraction decreases upon heating above 
Tf and disappears around 120 K. At high temperatures 
the main source of depolarization is the fast fluctuat- 
ing spins and the corresponding volume fraction, Af, in- 
creases linearly on warming (for 120 K< T < Tco)- This 
suggests that in this temperature range the short-range 
(CE-type) or dynamic CO dominates the muon depolar- 
ization, as found for x = 0.5 ^(| and Lai-^Ca^MnOa 
[l7|. The zero-field-cooled magnetization decreases al- 
most linearly with temperature for 120 K< T < Tco due 
to the formation of short-range antiferromagnetic corre- 
lations in the CO region (Fig. OJa)). 

In the limit of fast fluctuations one can estimate the 
correlation time, r, of magnetic fluctuations by assuming 
that the static width of the field distribution is of the or- 
der of the average magnetic field sensed by muons at the 
lowest temperature (7 A1 A_B A1 ~ 2nv^). We estimate short 
correlation times associated with fast fluctuations in the 
range Tf=10~ 12 -10~ 13 s for all samples, in good agree- 
ment with other values reported for bilayer manganitcs 
0, 0. For x=0.68, the additional slowly fluctuating 
component observed between 70 K and 120 K fluctuates 
with a correlation time r s =10 _9 -10 -11 s. The existence 
of slow and fast dynamics between Tf < T < 120 K for 
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x=0.68 further sustains the phase separation scenario in 
which charge-ordered and ferromagnetic (or A-type or C- 
type AFM) regions may coexist in the NLRO region, as 
found in Lai-^Ca^MnOs |20| . 
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cluster formation and /or magnetic polarons |22t |23| or 
stripe-like structures [24j. 

In conclusion, our /iSR study of electron-doped bilayer 
manganites yields local magnetic information and pro- 
vides evidence for magnetic order, phase separation and 
spin freezing. The results are summarized in the phase 
diagram shown in Fig. 4(a). Spin freezing is observed as 
a maximum in 1/T\ for samples with £=0.52-0.75 [l2j . 
The spin frozen state thus exists in the region of the 
phase diagram with LRO (x=0.52-0.62) together with a 
LRO state (AFM for £=0.52) and also in the region for 
which neutron diffraction failed to detect any LRO. We 
note that the spin-frozen region of the phase diagram 
overlaps closely with a region in which reentrant charge- 
ordering has been suggested to exist yjj, [l5j, |25j, |2fiJ . We 
can quantify the degree of local magnetic inhomogeneity 
at low temperatures using 1 / {2irv (Fig. 4(b)) which 
is largest for the x = 0.68 composition in the NLRO re- 
gion jg. 
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FIG. 4: (a) Phase diagram of La2-2xSri + 2 a; Mn207 (0.5 < x < 
1). The charge ordering temperatures were obtained from 
magnetization data .121 . Solid and dashed lines are guides 
to the eye. (b) Degree of magnetic inhomogeneity at 1.8 K 
expressed by 1 / \2nv M T2) . 

The freezing mechanism begins with the formation of 
magnetic clusters at high temperatures. Thermal disor- 
der opposes this tendency, so that with decreasing tem- 
perature the regions of correlated spins become larger 
and the spin system is subdivided into independent mag- 
netic clusters of different size with a distribution of block- 
ing temperatures. If a cluster is blocked for a time greater 
than the time window of the /iSR technique, it will ap- 
pear frozen or quasistatic and hence A r \ K decrease re- 
flects a progressive blocking of clusters by cooling down 
towards Tf , where all clusters are blocked. Spin- freezing 
phenomena are found in numerous s yste ms, ranging from 
spin-glasses, such as Cu-Mn alloys |2l| to lightly-doped 
cuprates, La2- x Sr K Cu04 Q. In our layered manganites 
we find that short-range charge-ordered fluctuations are 
dynamic at high temperatures (T < Tqo) and become 
static below Tf (or in the reentrant charge-ordered re- 
gion). We note that there has been evidence for spin 
freezing in manganese perovskites ascribed to magnetic 



[1 
[2 
[3 
[4 
[5 
[6 
[7; 
[8 
[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 

[23 
[24 
[25 
[26 



E. Dagotto et al, Physics Reports 244, 1 (2001). 

M. Uehara et al., Nature 399, 560 (1999). 

J. M. De Teresa et al, Nature 386, 256 (1997). 

C. Niedermayer et al, Phys. Rev. Lett. 80, 3843 (1998). 

C. D. Ling et al, Phys. Rev. B 62, 15096 (2000). 

J. E. Millburn et al, Chem. Comm. 15, 1389 (1999). 
R. H. Heffner et al, Phys. Rev. Lett 81, 1706 (1998). 
J. H. Brewer et al, Hyp. Int., 63, 177 (1990). 
1/Ti was also extracted over At = 0.1 — 25 /is. 

D. N. Argyriou et al, Phys. Rev. B 61, 15269 (2000). 
T. Chatterji et al, Phys. Rev. B 61, 570 (2000). 

A. I. Coldea, PhD Thesis, Oxford University (2001). 

S. T. Bramwell et al. J. Phys.: CM 5, 53 (1993). 

S. Rosenkranz et al, Physica B, 312, 763 (2002). 

J. Q. Li et al, Phys. Rev. B 64, 174413 (2001). 

D. Andreica et al, Physica B 289, 65 (2000). 

K. H. Kim et al, Phys. Rev. B 62, 11945 (2000). 

T. G. Perring et al, Phys. Rev. Lett. 78, 3197 (1997). 

J.-W. Feng et al, Phys. Rev. B 63, 180401 (2001). 

R. H. Heffner et al, Phys. Rev. B, 63, 094408 (2001). 

A. W. Hunt et al, Phys. Rev. Lett. 82, 4300 (1999). 

G. Papavassiliou et al, Phys. Rev. Lett. 87, 177204 

(2001). 

G. Allodi et al, Phys. Rev. Lett. 87 127206 (2001). 
T. Hotta, Phys. Rev. Lett. 86, 4922 (2001). 
Q. Yuan et al, Phys. Rev. Lett. 83, 3502 (1999). 
J. Dho et al, J. Phys.: Cond. Matt. 13, 3655 (2001). 



